//I/ s ¥ 

Cj4/< 


7 7£>3 2 ^ 

03 3 S 


THERMAL CONTROL SYSTEMS FOR 
LOW-TEMPERATURE HEAT REJECTION 
ON A LUNAR BASE 


K. R. Sridhar 
Principal Investigator 
and 

Matthias Gottmann 
Graduate Research Assistant 


Department of Aerospace and Mechanical Engineering 
The University of Arizona 
Tucson, AZ 85721 


Semiannual Status Report 
for Grant NAG5-1572 
from NASA Goddard Space Flight Center 


February 19 9 2 


(NASA-CR-190063) THERMAL CONTROL SYSTEMS 
FOR LOW-TEMPERATURE HEAT REJECTION ON A 
LUNAR BASE Semiannual Status Report 
(Arizona Univ.) 35 p CSCL 06K 

G3/54 


N9 2-20269 


Unci as 
0077032 


Nomenclature 


T 

— temperature 

L 

P 

— pressure 

A 

h 

— enthalpy 

H 

s 

— entropy 

S 

P 

— density 

(Jy 

m 

— mass flow rate 

d 

V 

— volumetric flow rate 

l 

Q 

— heat load 

f 

q 

— specific heat load 

V 

w 

— compressor power 

Be 

w 

— specific compressor power 

Cp 

COP 

— coefficient of performance 

p 

V 

— efficiency 

a 

m 

— specific mass 

e 

M 

— mass 



rejection loop pipe length 

radiator area 

radiator equivalent height 

radiator array distance to lunar base 

material strength 

inner pipe diameter 

pipe length 

friciton factor 

fluid velocity 

Reynolds number 

specific heat at constant pressure 

dynamic viscosity 

Stefan Boltzmann constant 

emissivity 
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1 Introduction 


One of the important issues in the lunar base architecture is the design of a Thermal Control 
System (TCS) to reject the low- temperature heat from the base. The TCS ensures that 
the base and all the components inside are maintained within the operating temperature 
range. The temperature of the lunar surface peaks to 400K during the 336 hour lunar day 
and heat rejection from the base under such conditions is a technically challenging task. 
Prior studies have shown that the overall mass of a TCS and its power supply for lunar 
base applications can be significant. The single largest fraction of the overall mission cost 
for any space mission is associated with the initial launch, which continues to be in the 
vicinity of $ 5000/lb from Earth to LEO. The reduction of lift mass at launch is a key driver 
in reducing the overall cost of space missions. In order to find the lowest mass for the TCS, 
several options have been proposed. One option would be to store the waste heat deep in 
lunar regolith [1], A piping system working as a heat exchanger has to be buried in the soil. 
The technical difficulties and uncertainties associated with large scale excavation on the 
Moon, and lack of knowledge about the thermal properties of lunar regolith are primary 
reasons for not pursuing this path presently. A significant portion of the total mass of 
the TCS is due to the radiator. In order to reduce the mass of the radiator, the concept 
of shaded light-weight radiators have been proposed [2], Shading the radiator from the 
sun and the hot lunar soil could decrease the radiator operating temperature significantly. 
This technology requires shades built of specular surfaces. The degradation of the radiator 
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Figure 1: Schematic Diagram of a TCS using a Heat Pump 

surface properties with age in a lunar environment is not known. At least for the initial 
cases, the prudent approach would be to employ systems that rely on proven technology. 
The concept of using a heat pump fits this bill. In this concept, energy in the form of 
heat or work, is supplied to the heat pump which collects heat from the low-temperature 
source (the lunar base) and delivers it at a higher temperature to the radiator. The mass 
of the radiator dissipating the high temperature heat would be significantly lower than 
one operating without a temperature lift. A simplified block diagram of this concept is 
illustrated in figure 1. 

Heat pumps have been in use for terrestrial applications for a long time. Refrigeration 
devices utilizing a thermodynamic cycle are essentially heat pumps. A vapor compression 
cycle involving two constant pressure and two adiabatic processes is the most widely used. 
It is also called a Rankine cycle and requires shaft work. Absorption cycles on the other 
hand are heat driven and do not require high quality shaft power. The Stirling cycle consists 
of two isothermal and two constant volume processes and promises a better efficiency than 
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the Rankine cycle. Theoretically, it reaches the same efficiency as the optimal Carnot cycle 
but the processes are technically difficult to realize. Today, Stirling cycle coolers are used 
in cryogenic applications. Experiments using this cycle for residential heat pumps show 
promising results [3], [4], but but these heat pumps are yet to become technically reliable. 

To optimize the mass of the heat pump augmented TCS, all promising options have 
to be evaluated and compared. During these preliminary comparison studies, considerable 
care is given to optimizing system operating parameters, working fluids and component 
masses. However, to keep this preliminary study simple and concise, the following aspects 
are not considered presently. The systems are modeled for full load operation and the im- 
plications and power penalties at off-design and partial load conditions are not considered. 
However, it is realized that the surface temperature of the lunar regolith varies consider- 
ably during the lunar day as shown in figure 2 . This -variation in the regolith temperature 
would indicate that the temperature lift and the load of the heat pump would vary as 
a function of the time of the day. For this reason, the performance of the heat pump 
at partial load conditions is important and will be studied in detail in the future. Also, 
redundancy requirements are not considered presently. While evaluating system mass, the 
control components are not accounted for since the variation in their masses for the various 
cycles and working fluids would not be large. Issues such as these will be studied in detail 
during the design of the actual system. 
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2 Rankine Cycle Heat Pump 


A study was conducted to optimize a heat pump operating on a Rankine cycle. The details 
are presented in this and the following sections 


2.1 Cooling Load 

To estimate the cooling load a closed system analysis was performed on the lunar base. 
Energies crossing the boundaries are electrical power supply, conduction through walls, and 
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heat removed by the acquisition loop. Within the system, heat generation can occur due 
to human metabolic activity. The electrical power input for a first stage base is estimated 
to be between 50 and 100 kW, more likely lOOkW ([1],[6] , [7]). The conduction through the 
walls depends on the insulation. Without significant mass penalties it is possible to reduce 
heat gains or losses to a very small fraction of the electrical input, hence it is neglected. 
Based on food consumption, a crew member produces an average of about 150W. For a 
crew of six to eight members, the total heat generation would again be negligible, compared 
to the electrical input. Therefore, the cooling load (the heat removed by the acquisition 
loop) can be equated to the electrical input into the base. Stated differently, this implies 
that all electrical input will finally be dissipated as heat. The value for the cooling load is 
fixed at lOOkW for this study. When further details about the design and activities of the 
base are known, these assumptions can be revisited and refined if necessary. 

2.2 The Acquisition Loop 

The acquisition loop collects the excess heat from the lunar base and transports it to the 
heat pump. The acquisition loop consists of cold plates and a network of connecting pipes. 
The heat is transported by a single phase fluid. Since the coolant in the acquisition loop 
circulates in the habitation module nontoxicity is a necessity from safety considerations. 
Water with certain trace additives to depress its freezing point would be a good candidate. 
For this study, it was decided that one cooling loop operating at a single pre-designed 
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temperature would be used. This temperature was chosen to be 275K (the lower of the 
two Space Station cooling loop temperatures). The variation in the temperature of the 
coolant has to be small enough to provide isothermal cooling for small variations in the 
load, yet large enough to keep the coolant flow rate within reasonable limits. The mass 
flow rate in the acquisition loop is rh = Q coo i/(cpZ^T). If water with trace additives were 
used as the coolant and the temperature variation in the acquisition loop were taken to be 
5K, the mass flow rate in the acquisition loop would be 4.8kg/s. 

2.3 The Heat Pump 

Two different heat pump setups were studied: In the first setup termed Case A the heat 
pump is directly connected to the rejection loop. In this case the condenser of the heat 
pump and the radiator are essentially the same device. The refrigerant circulating in the 
heat pump condenses and rejects heat in the radiator. An alternative setup is to decouple 
the heat pump and the rejection loop with a heat exchanger (Case B). Detailed analysis 
of the two cases and their pros and cons will be discussed in the following sections. 


2.4 Heat Pump Connected Directly to the Rejection Loop (Case A) 


A simplified schematic of a heat pump configuration directly connected to the rejection 
loop is illustrated in figure 3 The main parameters of interest in the design of a heat 
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Figure 3: Schematic Diagram of a Heat Pump directly connected to the Rejection Loop 


pump used for cooling are the input heat load (Q coo j) and its temperature (Tj ow ), the 
temperature lift, and the coefficient of performance (COP). The COP of a heat pump is 


defined as 


COP = 


Qcool 

W 


where W is the power consumed by the heat pump. 


2.4.1 The Compressor 


Figure 4 illustrates the Rankine cycle process in a pressure-enthalpy diagram. The working 
fluid in the vapor state is compressed from P\ to P 2 . Ideally this process would be isentropic 
(2S). Due to irreversibilities, the process is nonisentropic 


^ideal 


his — hi 






^real = h 2 - h x 

— ^ 2S ~ hi 
^comp ~ h2 _ hl 


where: P is the pressure, h is the specific enthalpy and the subscripts refer to the states 
in figure 4. It is assumed that the compression would be performed in a single stage, 
in order to limit the number of free parameters. Customarily airplane cooling systems 
utilize multistage compression [8], but there is no intercooling between the stages. Hence, 
effectively the compression can be modeled as a single stage process. The properties of 
the refrigerant used for the calculation are obtained from [9] and a FORTRAN77 code 
developed in-house [10]. Deviation from the ideal behavior in the compression occur due 
to mechanical, electrical (motor), and electronic (controller) inefficiencies, and fluid friction. 
The values for the efficiencies in state-of-the-art aircraft cooling equipment were obtained 
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from [11], and are as follows: 7/ mech = 0.95, Electrical = 0.94, Electronic = 0-91, and 
7/£ u jd = 0.75. The overall efficiency of such a compressor would be the product of all 
four efficiencies, about 61 percent in this case. The excess energy supplied to overcome 
these inefficiencies will be converted to heat. Since the compressor would operate in a 
high vacuum environment, radiation to the environment and convection of the heat by 
the vapor flow inside are the only heat rejection mechanisms. An upper limit for the 
amount of heat rejected by radiation can be obtained by modeling the compressor as a 
black cube, 0.25m side, at 400K. Both, the surface properties and the temperature of the 
surface are deliberately chosen to be well above the actual values to obtain a conservative 
estimate . This heat is small compared to the compressor input power and can be neglected. 
Therefore, it can be assumed that all the energy supplied to the compressor will be used 
to compress and heat the refrigerant. It should be noted that the temperature of the 
compressor can be maintained within operating limits by the use of a cold plate. However 
this would not be required since the working fluid could convectively remove the excess 
heat from the compressor. 

The next step is a mass estimate for the compressor. In aircraft cooling, the compressor 
mass is assumed to be proportional to the cooling load. One pound (0.454kg) per kilowatt 
is the value suggested [11]. Our mass optimization is performed at maximum cooling load 
(lOOkW) and this value remains unchanged in the optimization. The heat pump output 
temperature and hence the total heat rejected by the heat pump is varied. Since the 
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assumption of compressor mass being proportional to the cooling load would lead to an 
unrealistic constant mass estimate in our case, it was modified as follows. A proportionality 
was assumed between compressor mass and the heat pump output heat, which is the sum 
of the input heat and compressor power. The proportionality constant was arrived at as 
follows. The reference temperatures to obtain the proportionality factor, Tj^gh = 380K 
Tjow = 275K, are values typical for an aircraft cooling system. For these temperatures 
and R717 as refrigerant, the heat pump overall COP is 0.805. With this value we get: 


^comp — 


_ Afpomp __ 


^reject 


^cool m comp] 

COP 


= 0.202 


k g 

kW 


where m C omp is the compressor mass in kg per kW rejected heat, A/comp the actual 
compressor mass in kg, m^omp the compressor mass in kg per kW cooling load. 


2.4.2 Discharge and Return lines to and from the Radiator 

At state 2 in figure 4 the refrigerant is in the superheated state. The length of the discharge 
line depends on the layout of the lunar base and how the radiators are spatially configured. 
The discharge line has to connect all the radiators to the compressor. Assuming the 
radiators are of constant height, it is reasonable to take the pipe length to be proportional 
to the radiator area, i.e. L — S + A/H, where L is the length that will be used to determine 
the pressure drop, A the radiator area, H is an “equivalent height” of the radiator, and S 
is the distance from the lunar base to the radiator array. Figure 5 depicts schematically 
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Figure 5: Schematic of Radiators and Rejection Loop Piping 


the setup of the radiators and the piping. The term “equivalent height” is used because 
this dimension can differ from the actual height to accommodate bends in the piping or a 
spacing in between the radiators. The pipe length was defined as the distance from the heat 
pump to all radiators. The complete rejection loop is 2 L. The pressure drop in the piping is 
a function of the pipe diameter, and is determined based on handbook recommendations for 
good design practice [12]. The pressure drop in the discharge line, the radiator (condenser) 
and the return line is taken to be the equivalent of a IK temperature drop. It is important 
for the thermodynamic model that this pressure drop be small enough in order for it to not 
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affect the overall efficiency. The fixing of the total pressure drop also allows the designer 


to decouple the pipe sizing from the thermodynamic evaluation of the heat pump. The 
pressure drop is split such that one half of it occurs in the condenser and the rest is in the 
discharge and return lines. The friction losses in discharge and return lines are determined 
based on the optimization of the pipe masses. The frictional pressure drop, (A P)j = 




-2 


d is the pipe diameter, 


where the friction factor for smooth pipes / = 

/ is the length of the pipe, v the fluid velocity, p the fluid density, and Be the pipe Reynolds 
number. The total mass is the sum of the mass of the pipe and the mass of the fluid in 
the pipe. The tube thickness is computed based on a factor of safety of three. A minimum 
thickness of 0.5mm is also required. The density of the piping material is based on a light 
weight, high strength aluminum alloy. Should such an alloy be chemically incompatible 
with the refrigerant of choice, the inside of the pipes can be surface coated to take care of 
the problem. The masses are: 


M, 


pipe 


^Ppipe^ 

2<T y,pipe 


■^fluid 


4 


where <7y pjp e is the allowable (design) stress for the pipe material. 

Between states 2 and 3 the superheated vapor is cooled in the radiator. Ideally this 
process can be modeled as an isobaric process, but due to pipe friction a small pressure 
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drop would occur. Between states 3 and 4, the refrigerant is condensed to saturated liquid. 
A finite pressure drop occurs in the condenser. The mass estimate for the condenser will be 
discussed in the radiator section. The heat to be rejected by the radiator, <? re j ec t = h^ — h^. 
From state 4 the saturated liquid is sent from the radiator to the throttle valve located at 
the evaporator inlet, through the return line. The sizing of the return line is based on the 
same guidelines described for the discharge line. 

2.4.3 Evaporator and Throttle Valve 

Between states 4 and 5 the fluid is adiabatically throttled. The mass of the throttle valve 
is negligible compared to the mass of the other components of the heat pump. Between 
states 5 and 1, the refrigerant absorbs heat from the primary coolant circulating in the 
lunar base. The temperature difference between the primary coolant in the acquisition loop 
and the boiling refrigerant in the evaporator is 5K. The heat removed, <7 coo ] = hi — h 5 . 
The mass of the evaporator is obtained based on a suggested value of 2.72kg/kW [1]. 

2.4.4 Refrigerant 

One of the important issues is the choice of refrigerant as the working medium for the 
Rankine cycle. The refrigerants that are commonly used in terrestrial and aerospace ap- 
plications, Rll, R12, R113, R114, and R717, were considered [8]. R113 and R114 were 
eliminated from the list of potential refrigerants due to the possibility of condensation of 
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Figure 6: Rankine Cycle for R114. 7] ow = 270K, = 360K 
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the vapor in the compressor (figure 6). Such condensation would be detrimental to the life 
of the compressor. Among the remaining refrigerants, the selection was narrowed down 
to Rll and R717, because R12 has a lower COP and a lower critical temperature (R717: 
T^rit = 407K, Rll: T cr ^ = 474K, R12: T cr j{. = 385K). The P — h diagrams for R717 and 
Rll are shown in figures 4 and 7 respectively. Safety considerations give an edge to Rll 
due to its nontoxicity and noninflamability, but R717 offers lower overall system mass. The 
thermodynamic properties of the refrigerants were obtained using the analytical functions 
suggested in [9]. The COP can be expressed in terms of the specific enthalpies and for the 


Rankine cycle: 


COP = 
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Enthalpy [kJ/kg] 

Figure 7: Rankine Cycle for Rll. Tj ow = 270K, T^igh = 360K 


Table 1 illustrates the COP calculation for a condenser temperature of = 380K. The 

overall COP was computed as a function of the condenser temperature and is plotted in 
figure 8. 

2.4.5 Implementation of Heat Pump and Piping Model 

Values for COP and the mass of the piping were computed and tabulated for varying 
rejection temperature using the models presented above. These tabulated values were 
imported into the spreadsheet and linearly interpolated where necessary. 
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Refrigerant: R717 


state 

T 

P 

h 

s 

P 


[K] 

[MPa] 

[kj/kg] 

[k J /kg] 

[kg/m 3 ] 

1 

270 

0.381 

1584 

6.046 

3.088 

2 

626 

7.27 

2419 

6.615 

24.89 

3 

380 

7.14 

1541 

4.788 

67.2 

4 

380 

7.14 

893 

3.080 

436.5 

5 

270 

0.381 

893 

3.483 

6.725 

COP = 0.829 


Refrigerant: Rll 


state 

T 

P 

h 

s 

P 


[K] 

[MPa] 

[kJ/kg] 

[kj/kg] 

[kg/m 3 ] 

1 

270 

0.035 

249 

0.9581 

2.18 

2 

448 

0.964 

350 

1.051 

39.2 

3 

380 

0.945 

301 

0.9343 

49.7 

4 

380 

0.945 

158 

0.5565 

1255 

5 

270 

0.035 

158 

0.6180 

4.22 

COP = 0.914 


Table 1: Properties in the Refrigeration Cycle for = 380K 

2.5 Heat Pump Decoupled from the Rejection Loop with a 
Heat Exchanger (Case B) 


Connecting the heat pump directly to the radiator has inherent disadvantages. If the 
refrigerant used in the Rankine cycle is not suitable for a heat transport loop, it can be 
advantageous to separate the rejection loop from the heat pump with a heat exchanger. 
This configuration of a heat pump augmented TCS is shown in figure 9. From a system 
design perspective it is desirable to decouple subsystems that carry out different tasks. 
The decoupled case would provide for better and simpler control of the TCS during partial 
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Figure 9: Schematic Diagram of a Heat Pump decoupled from the Rejection Loop 


load conditions. On the other hand, a heat exchanger between the two loops will cause a 
temperature drop between the heat pump and the rejection loop and an associated mass 
penalty. To compensate for the temperature drop the heat pump has to deliver the output 
heat at a higher temperature and therefore operate at a lower COP. If the same fluid were 
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used in the Rankine cycle and in the rejection loop, the only foreseeable advantage of the 
decoupled system would be the possibility of better and simpler control. However, other 
advantages could emerge if two different fluids were used. 

The thermodynamic and mass models for the heat pump with an output heat exchanger 
(Case B) differs only in few aspects from the models presented for Case A (section 2.4). 
Only these differences will be discussed in this section. 

2.5.1 Condenser 

In case B the condenser is a heat exchanger that decouples the rejection loop from the 
heat pump. Both fluids undergo phase change in this heat exchanger. For a mass estimate 
the value quoted in [1], 2.72 kg/kW was used. The thermodynamic performance of the 
condenser is characterized by a pressure drop in each loop (heat pump and rejection loop) 
and a temperature difference between both sides. Similar to the acquisition side, the 
temperature difference is set to 5K. Consistent with case A, the pressure drop has to be 
small enough, not to affect the heat pump performance. A pressure drop equivalent to IK 
temperature drop has been assigned to the condenser. 

2.5.2 Rankine Cycle Analysis 

The cycle evaluation follows the same path outlined for Case A. The efficiencies and pres- 
sure drops of the heat pump components are also the same as in Case A. The COP as a 
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Figure 10: Comparison of COP from Cycle Analysis and Approximation for Rll and R717. 
^evaporator = 270K 


function of the output temperature T^g^ was computed with a FORTRAN77 program 
using real fluid properties from [9]. The implementation of this COP(T ) into the spread- 
sheet was realized with an approximate analytical function. For each refrigerant a fourth 
order polynomial was fitted to the data computed with the FORTRAN77 code. The re- 
sulting approximation yields an error of less then 0.3 percent for output temperatures from 
T^igh = 320K to T^jg^ = 390K. Figures 10a and 10b show a comparison between the real 
fluid model and the polynomial approximation. It can be seen that the results are almost 
identical. 


2.5.3 Rejection Loop 


The decoupled rejection loop would require a pump to circulate the coolant fluid. This 
pump and the power penalty associated with it have to be incorporated in the mass estimate 
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and optimization. The pump mass is estimated using a formula quoted by Dexter and 


Haskin [13] 


Afpump 



where: m is the mass flow rate in lb/hr, p is the density of the fluid in lb/ft 3 . The power 
required for a liquid pump can be readily computed from 


Wpump — 


APV 

7/pump 


where A P is the pressure differential across the pump, V the volume flow rate, and ?/purnp 
the pump efficiency. A conservative ?7pump = 0.25 as suggested by [13] was used. The 
pressure drop was determined with the formulas presented for case A. The pipe thickness is 
again determined based on the hoop stress or 0.5mm, whichever is larger. Masses included 
in the estimate are due to pipes, coolant, pump, and the power supply. The decoupled 
rejection loop does not affect the heat pump COP. The minimum mass for the loop may be 
achieved by balancing pipe mass and the power penalty. This approach results in optimum 
mass when the pipe diameters are relatively small and the pressure drop is large. However, 
a large pressure drop in the vapor line would result in a large temperature drop and this is 
accompanied by an increase in the radiator area and mass. While the pressure drop in the 
liquid line can be compensated with the pump, if the pressure drop gets large, the pumping 
power will become significant and add to the total heat rejection load.. Therefore, the mass 
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estimate for the piping has to be computed based on a limited pressure drop. Here again, 
the pressure drop is specified in terms of an equivalent temperature drop and is set to 0.5K 
in the vapor and 1.0K in the liquid line. These values are chosen based on recommended 
design practice [12]. The cooling fluid of choice is ammonia, which already demonstrated 
its good performance as a heat transport fluid in Case A. The toxicity of ammonia will not 
be a concern for the rejection loop as it is outside the habitation modules. 

In Case A the piping mass was determined together with the heat pump estimate be- 
cause they are coupled. Assuming values for the radiator height and distance from the base, 
Case A yielded a model where the piping mass depends solely on the rejection temperature. 
For Case B a model that makes use of the decoupling of heat pump characteristics and 
the rejection loop was sought. For a given refrigerant and specified pressure drops in the 
liquid and vapor lines, the rejection loop mass depends on three parameters: rejection heat 
load Q re j ec t) rejection temperature T re j ec ^, and pipe length L. Using the thermodynamic 
properties from [9] the mass model was implemented into a FORTRAN77 code. Figures 
11 and 12 show results obtained with the code. For use with a spreadsheet software it is 
desirable to obtain an analytical expression for the mass. This was realized with a polyno- 
mial which is second order in temperature, second order in length and linear in rejection 
heat load: 

"piping = £X:i>*rW 
i~0 j~0 k = 0 

The coefficients were determined with a least square error fit. The approximation is valid 
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Figure 11: Mass of the Liquid Piping 


in the following range: 340K < T re j ec t < 380K, 150kW < Q re j ec t < 250k W, and 100m < 


L < 400m. The maximum error of the approximation is three percent. 


2.6 Radiator Considerations 

The function of the radiator is to reject the waste heat from the base. The heat rejected 
by the radiator is given by Q — ^ e7 7 cr (^r e j ec t — ^sink) w ^ ere e * s emissivity, 77 the 
fin efficiency, T re j ec t and the radiator and sink temperatures. The estimated sink 

temperature for a vertically mounted radiator at the lunar base is 32 IK [6]. Most reviewed 
sources suggest e = 0.8, and rj = 0.7. Several estimates for the mass of a radiator are avail- 
able in the literature ([1], [6], [14], [15]). The mass of a radiator is taken to be proportional 
to its area and recent publications recommend a value of 5kg/m 2 for a one sided radiator. 
The vertical radiator is two sided and hence a mass estimate of 2.5kg/m 2 is assumed. The 
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Figure 12: Mass of the Vapor Piping 


heat to be rejected is the cooling load of the base plus the power consumed to operate the 
heat pump. 


2.7 Power Supply 


The heat pump consumes power to achieve the desired temperature lift. The capacity of the 
lunar base power station needs to be greater than otherwise, to account for this additional 
power consumption. It is reasonable to assume that the additional mass penalty would be 
proportional to the power supplied to the heat pump. A review of the literature shows 
that there is no consensus on the mass penalty ([1], [6], [16], [17]). The values quoted lately 
are in the neighborhood of 30kg/kW for nuclear units of the SP100 class [15]. This value 
will be used in our studies. 
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3 Results 


The overall mass optimization was performed in a spreadsheet. The heat pump output 
temperature lift and hence the radiator temperature was varied and the variation of the 
masses of the components and the TCS were computed using the mass models described 
here. For the coupled TCS configuration, Case A, the analysis were performed for two 
working fluids, Rll and R717 and the overall TCS mass variation as a function of radiator 
temperature is shown in figure 13a. Similar analysis were performed for a the decoupled 
configuration, Case B, and the results are shown in figure 13b. For Case B, Rll and R717 
were used as working fluids for the heat pump, but R717 was used in the rejection loop 
due to to its superior heat transport characteristics. 

When Rll is used as the working fluid for the heat pump, the optimal TCS mass is 
6108kg at a radiator temperature of 371K for the coupled situation, Case A. For Case B, 
the optimal TCS mass is 5940kg at a radiator temperature of 362K. The radiator mass 
in Case B is higher than for Case A due to its lower operating temperature. Also, the 
presence of the heat exchanger between the heat pump and the rejection loop adds extra 
mass to the Case B scenario. In spite of these mass penalties the optimal TCS system 
mass for Case B is lower than that for Case A. This is due to the huge reduction in the 
rejection loop piping mass for Case B. When R717 is used as the working fluid in the heat 
pump, the optimal mass of the TCS is 5515kg at a radiator temperature of 362K for Case 
A. For Case B the corresponding values are 6392kg and 360K respectively. It is obvious 
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that Case B is more massive than Case A since the radiator temperature for Case B is 
lower and also it has an additional heat exchanger. 

The mass of the individual components for Cases A and B are shown graphically in 
figures 14 and 15 for a range of radiator temperatures and is listed in tables 2-5 for the 
optimal radiator temperatures. 

Among the cases considered R717 coupled TCS configuration offers the least mass of 
5515kg. The best decoupled configuration would involve Rll as the working fluid for the 
heat pump and R717 as the working fluid for the rejection loop. The optimal mass for this 
configuration, as stated earlier, is 5940kg. In spite of the additional mass, the decoupled 
system may offer the advantage of simpler and more reliable system control during partial 
load operation. This issue will be examined in detail in the future. 


b ° b ? 



Figure 13: Overall TCS Mass as a Function of Tj-gjg^ 
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ACQUISITION LOOP 



Cooling Load 

Q cool 

100 

kW 

Cooling Temperature 

Tcool 

275 

K 

HEAT PUMP 



Temperature Drop HXi n 

A Tuxin 

5 

K 

Input Temperature 

Tlow 

270 

K 

Output Temperature 

Thigh 

362 

K 

Heat Pump Efficiency 

COP 

1.11 


Compressor Power 

W 

90.2 

kW 

Rejection Heat Load 

Q reject 

190.2 

kW 

Evaporator Spec. Mass 

Wlevap " 

2.72 

kg/kW 

Compressor Spec. Mass 

'Wlcomp 

0.202 

kg/kW 

Evaporator Mass 

evap 

272 

kg 

Compressor Mass 

M CO mp 

18.2 

kg 

Heat Pump Mass 

MjqP 

290.2 

kg 

POWER SUPPLY 



Power Supply Spec. Mass 

'H'lpower 

30 

kg/kW 

Power Penalty 

JMpower 

2707 

kg 

REJECTION LOOP 



Pipe Mass 

Mpip e 

278.2 

kg 

RADIATOR 



Rejection Temperature 

T r eject 

362 

K 

Sink Temperature 

Tsi n k 

320 

K 

Fin Efficiency 

V 

0.7 


Emissivity 

e 

0.8 


Radiator Area 

A 

895.9 

m 2 

Radiator Spec. Mass 

TYlrad. 

2.5 

kg/m 2 

Radiator Mass 

Mrad 2239.8 

kg 

SYSTEM 



System Mass 

Mtcs 

5515 

kg 


Table 2: Optimum Component and TCS Mass for Case A with R717 
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ACQUISITION LOOP 


Cooling Load 

Qcool 

100 

kW 

Cooling Temperature 

Tcool 

275 

K 

HEAT PUMP 



Temperature Drop HX; n 

kTnxin 

5 

K 

Input Temperature 

T low 

270 

K 

Output Temperature 

Thigh 

371 

K 

Heat Pump Efficiency 

COP 

1.06 


Compressor Power 

w 

94.5 

kW 

Rejection Heat Load 

Q reject 

194.5 

kW 

Evaporator Spec. Mass 

m evap 

2.72 

kg/kW 

Compressor Spec. Mass 

comp 

0.202 

kg/kW 

Evaporator Mass 

-M eva p 

272 

kg 

Compressor Mass 

M c omp 

19.1 

kg 

Heat Pump Mass 

Mhp 

291.1 

kg 

POWER SUPPLY 



Power Supply Spec. Mass 

power 

30 

kg/kW 

Power Penalty 

Mpower 

2836 

kg 

REJECTION LOOP 



Pipe Mass 

M pipe 

1170.5 

kg 

RADIATOR 



Rejection Temperature 

T r eject 

371 

K 

Sink Temperature 

Tsink 

320 

K 

Fin Efficiency 

7 

0.7 


Emissivity 

e 

0.8 


Radiator Area 

A 

724.2 

m 2 

Radiator Spec. Mass 

TTlrad 

2.5 

kg/m 2 

Radiator Mass 

Atrad 

1810.5 

kg 

SYSTEM 



System Mass 

Mtcs 

6108 

kg 


Table 3: Optimum Component and TCS Mass for Case A with Rll 
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ACQUISITION LOOP 



Cooling Load 

Qcool 

100 

kW 

Cooling Temperature 

Tcool 

275 

K 

HEAT PUMP 



Temperature Drop HX; n 

A Tux in 

5 

K 

Temperature Drop HX out 

HXout 

5 

K 

Input Temperature 

Tlow 

270 

K 

Output Temperature 

Thigh 

365 

K 

Heat Pump Efficiency 

COP 

1.06 


Compressor Power 

w 

94.0 

kW 

Rejection Heat Load 

Q reject 

194.0 

kW 

Evaporator Spec. Mass 

TTlevap 

2.72 

kg/kW 

Condenser/HX Spec. Mass 

771 condL 

2.72 

kg/kW 

Compressor Spec. Mass 

771 comp 

0.202 

kg/kW 

Evaporator Mass 

M eV ap 

272 

kg 

Condenser/HX Mass 

Mcond 

527.8 

kg 

Compressor Mass 

Mcomp 

19.0 

kg 

Heat Pump Mass 

Mhp 

818.8 

kg 

POWER SUPPLY 



Power Supply Spec. Mass 

TPpower 

30 

kg/kW 

Power Penalty 

power 

2821.4 

kg 

REJECTION LOOP 



Liquid Pipe Mass 

liquid 

213.3 

kg 

Vapor Pipe Mass 

Myapor 

117.5 

kg 

Pipe Mass 

Mpip e 

330.8 

kg 

RADIATOR 



Rejection Temperature 

Treject 

360 

K 

Sink Temperature 

T s i n k 

320 

I< 

Fin Efficiency 

V 

0.7 


Emissivity 

e 

0.8 


Radiator Area 

A 

968.5 

m 2 

Radiator Spec. Mass 

TTlfad 

2.5 

kg/m 2 

Radiator Mass 

M ra d 

2421.3 

kg 

SYSTEM 



System Mass 

Mjcs 

6392 

kg 


Table 4: Optimum Component and TCS Mass for Case B with R717 
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ACQUISITION LOOP 



Cooling Load 

Qcool 

100 

kW 

Cooling Temperature 

T C ool 

275 

K 

HEAT PUMP 



Temperature Drop HX; n 

ATffXin 

5 

K 

Temperature Drop HX ou t . 

^Tnxout 

5 

K 

Input Temperature 

Ti ow 

270 

K 

Output Temperature 

Thigh 

367 

K 

Heat Pump Efficiency 

COP 

1.14 


Compressor Power 

w 

87.7 

kW 

Rejection Heat Load 

Q reject 

187.7 

kW 

Evaporator Spec. Mass 

TYlevap 

2.72 

kg/kW 

Condenser/HX Spec. Mass 

Wlcond 

2.72 

kg/kW 

Compressor Spec. Mass 

'Tft'comp 

0.202 

kg/kW 

Evaporator Mass 

Mevap 

272 

kg 

Condenser/HX Mass 

Mond 

510.6 

kg 

Compressor Mass 

M C omp 

17.7 

kg 

Heat Pump Mass 

Mhp 

800.3 

kg 

POWER SUPPLY 



Power Supply Spec. Mass 

m power 

30 

kg/kW 

Power Penalty 

M power 

2631.3 

kg 

REJECTION LOOP 



Liquid Pipe Mass 

M liquid 

193.5 

kg 

Vapor Pipe Mass 

M V apor 

104.8 

kg 

Pipe Mass 

Mpip e 

298.3 

kg 

RADIATOR 



Rejection Temperature 

T-r eject 

362 

K 

Sink Temperature 

T s ink 

320 

K 

Fin Efficiency 

V 

0.7 


Emissivity 

e 

O.S 


Radiator Area 

A 

884.2 

m 2 

Radiator Spec. Mass 

T^rad 

2.5 

kg/m 2 

Radiator Mass 

M r ad 

2210.4 

kg 

SYSTEM 



System Mass 

MjCS 

5940 

kg 


Table 5: Optimum Component and TCS Mass for Case B with Rll 
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